OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Phosphoproteomic Profiling of Human Myocardial ^ 
Tissues Distinguishes Ischemic from Non-lschemic End cros^Mark 
Stage Heart Failure 

Matthew A. Schechter^ ^ Michael K. H. Hsieh^ ^ Linda W. Njoroge\ J. Will Thompson^ Erik J. Soderblom^ 
Bryan J. Feger\ Constantine D. Troupes\ Kathleen A. Hershberger^'^, Olga R. Ilkayeva^ 
Whitney L. Nagel\ Gina P. Landinez\ Kishan M. Shah\ Virginia A. Burns^, Lucia Santacruz\ 
Matthew D. Hirschey^'"*, Matthew W. Foster^ Carmelo A. Milano\ M. Arthur Moseley^ 
Valentino Piacentino lll\ Dawn E. Bowles^* 

1 Department of Surgery, Duke University Medical Center, Durham, North Carolina, United States of America, 2 Duke Proteomics Core, Duke University Medical Center, 
Durham, North Carolina, United States of America, 3 Sarah W. Stedman Nutrition and Metabolism Center, Duke University Medical Center, Durham, North Carolina, United 
States of America, 4 Department of Medicine, Duke University Medical Center, Durham, North Carolina, United States of America, 5 Duke Translational Research Institute, 
Duke University Medical Center, Durham, North Carolina, United States of America, 6 Division of Pulmonary, Allergy and Critical Care, Medicine, Duke University Medical 
Center, Durham, North Carolina, United States of America 



Abstract 

The molecular differences between ischemic (IF) and non-ischemic (NIF) heart failure are poorly defined. A better 
understanding of the molecular differences between these two heart failure etiologies may lead to the development of 
more effective heart failure therapeutics. In this study extensive proteomic and phosphoproteomic profiles of myocardial 
tissue from patients diagnosed with IF or NIF were assembled and compared. Proteins extracted from left ventricular 
sections were proteolyzed and phosphopeptides were enriched using titanium dioxide resin. Gel- and label-free nanoscale 
capillary liquid chromatography coupled to high resolution accuracy mass tandem mass spectrometry allowed for the 
quantification of 4,436 peptides (corresponding to 450 proteins) and 823 phosphopeptides (corresponding to 400 proteins) 
from the unenriched and phospho-enriched fractions, respectively. Protein abundance did not distinguish NIF from IF. In 
contrast, 37 peptides (corresponding to 26 proteins) exhibited a >2-fold alteration in phosphorylation state (p<0.05) when 
comparing IF and NIF. The degree of protein phosphorylation at these 37 sites was specifically dependent upon the heart 
failure etiology examined. Proteins exhibiting phosphorylation alterations were grouped into functional categories: 
transcriptional activation/RNA processing; cytoskeleton structure/function; molecular chaperones; cell adhesion/signaling; 
apoptosis; and energetic/metabolism. Phosphoproteomic analysis demonstrated profound post-translational differences 
in proteins that are involved in multiple cellular processes between different heart failure phenotypes. Understanding the 
roles these phosphorylation alterations play in the development of NIF and IF has the potential to generate etiology-specific 
heart failure therapeutics, which could be more effective than current therapeutics in addressing the growing concern of 
heart failure. 
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Introduction 

Despite earlier diagnosis and improved therapy, heart failure 
(HF) continues to be a major health concern. In 2012, 5.7 million 
Americans were diagnosed with HF [1]. The lifetime risk of 
developing HF after age 40 is 20%, with the annual incidence 
approaching 10 per 1000 people after age 65 [2] . More than half 
of all HF patients will die within a 5-year period of being 



diagnosed [3,4]. These statistics are complicated by the fact that 
HF is a complex, multi-faceted disease that presents in two major 
forms: ischemic and non-ischemic. 

Ischemic HF describes significantly impaired left ventricular 
function resulting from reduced blood supply to the heart muscle. 
In contrast, the reduced left ventricular function seen in non- 
ischemic HF has a range of etiologies, including congenital, 
infectious agents, autoimmune, and idiopathic. However, the 
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current standard of care treats all cases of HF similarly, regardless 
of etiology. The treatment options for advanced heart failure are 
limited to implantation of a ventricular assist device to mechan- 
ically unload the heart, heart transplantation, or palliation with 
continuous intravenous inotropic support. These options, howev- 
er, are associated with high morbidity and mortality, highhghting 
the continued need for improving HF treatment. 

One such avenue of exploration is etiology-specific treatment. 
Such precise therapy would require an enhanced understanding of 
the molecular differences between heart failure phenotypes. 
Uncovering these molecular differences in a systematic and 
comprehensive way is made possible by utilizing high throughput 
'omics profiling [5,6]. Genomics, proteomics, metabolomics, etc. 
have primarily focused on comparing the difiFerences between one 
HF etiology' and non-failing samples (reviewed by Gonzalez et al 
[7]). Recentiy, using two dimensional gel-based proteomics on 
tissue from diseased human hearts, Klawitter el al were able to 
demonstrate differences in the relative abundance of proteins in 
signaling pathways between ischemic or idiopathic cardiomyop- 
athies [8] . However, differences in protein quantity represent only 
one aspect of the molecular picture; many of the key proteins 
involved in signal transduction pathways are highly regulated by 
post translational modification, such as phosphorylation. Whether 
protein phosphorylation patterns differ between heart failure 
etiologies is unknown. Such knowledge may enable the develop- 
ment of better and more precise heart failure therapeutics and 
diagnostics. 

In the current im-cstigation, cardiac, tissues from a well- 
characterized and weU-preserved human heart tissue bank were 
subjected to titanium dioxide resin to enrich for phosphopeptides, 
and then analyzed by a bottom-up LC/MS/MS global proteomics 
approach. This approach revealed amino acid residue-specific 
phosphorylation patterns on 400 cardiac proteins, which were 
compared between the IF and NIF etiologies. This revealed, for 
the first time, cardiac disease-specific phosphorylation pattern 
variations on key proteins involved in various aspects of cardiac 
physiology. Understanding how distinct protein phosphorylation 
patterns impact specific heart failure etiologies will support the 
development of therapeutics that better treat heart failure. 

Materials and Methods 

Human Cardiac Tissue Acquisition and Tissue Repository 

The study protocol was approved by the Duke University 
Health System Institutional Review Board (DUHS IRB). The 
DUHS IRB, FWA #00009025, is duly constituted and complies 
with all U.S. regulatory requirements related to the protection of 
human research participants and the Guidelines of the Interna- 
tional Conference on Harmonization (ICH) as adopted by the 
U.S. Food and Drug Administration. Tissue samples used for this 
study were procured from the Duke Human Heart Repository 
(DHHR), which is a DUHS IRB appro\cri tissue repository. 
Samples were procured by the DHHR in accordance under an 
approved DUHS IRB protocol using written informed consent or 
a waiver of consent for discarded tissues. Some samples procured 
had been anonymized with all HIPAA identifiers removed, others 
were de-identified. Only DHHR personnel were pri\y to a key 
that could link the samples to patient information. No HIPAA 
information was provided with any of the samples used in this 
study. Human myocardium was acquired from the left ventricular 
free wall of explanted IF or NIF hearts following cardiac 
transplantation. Non-failing (NF) left ventricular tissue was 
acquired from donors whose hearts were not utilized for 
transplant, thus becoming available for research. After explanta- 



tion, transmural tissue samples were processed and stored as 
described in the Methods SI. 

Sample Preparation for Mass Spectrometry 

Heart tissue samples were homogenized and subjected to phase 
separation by addition of chloroform. Protein was precipitated 
from the organic layer, washed, sonicated, recovered by centrifu- 
gation, and re-suspended in mass spectrometry-compatible deter- 
gent (RapiGc-st, AVaters Corp., Milford, MA). Cardiac tissue 
homogenati's used for mass spectrometry were subjected to 
Bradford Assay for protein quantification. A 625 |J,g aliquot of 
protein (per sample) was subjected to reduction, alkylation, 
followed by overnight proteolysis with sequencing grade trypsin 
(Promega, Madison, WI). A 25 |J,g aliquot from each sample was 
used for unenriched proteomic analysis of protein expression in 
the heart tissue. This 25 |Xg was spiked with 1.25 pmol 
ADH1_YEAST digest (Massprep standard. Waters Corporation) 
as a surrogate standard prior to analysis. 

The remaining 600 |j,g of protein were then enriched for 
phosphopeptides using in-house packed Ti02 spin columns as 
previously described [9]. 

LC/MS/MS Data collection 

The sample cohort was randomized prior to LC/MS/MS 
analysis. Peptide digests obtained from each of the samples were 
analyzed in a label-free quantitative fashion using a nanoAcquity 
UPLC system coupled to a Synapt HDMS mass spectrometer 
(Waters Corp, Milford, MA) for unenriched peptide analyses and 
an LTC) Orbitrap XL (Thermo Fisher Scientific, Waltham, MA) 
for phosphopeptide analyses. 

LC-MS Data Processing 

Robust peak detection and label-free alignment of individual 
peptides across all sample injections was performed using the 
commercial package Rosetta Elucidator v3.3 (Rosetta Biosoftware, 
Inc., Scattk;, WA) with PeakTeller algorithm [10]. Details are 
described in supporting methods. The raw data in the form of a 
flat file (excel) Excel File SI is provided via this fink: https:// 
discovery.genome.duke.edu/ express/resources/3772/Schechter_ 
PLOSone_Supplemental_Data June20 1 4.xlsx. 

Statistical Analysis 

Basic statistical analysis was performed on both the unenriched 
(protein-level) and phosphopeptide (peptide-level) datasets in order 
to obtain candidate (phospho) proteins that were differentially 
expressed. Fold-changes were calculated for each failing group 
versus non-failing control, as the ratio of the average intensity 
between the groups; directionaUty of the ratio was established that 
positive fold-changes mean up-regulated in failing versus nonfail- 
ing control, and negative fold-changes mean down-regulated in 
failing versus nonfailing. P-valucs were calculated using an error- 
weighted ANOVA with Benjamini-Hochberg FDR correction for 
multiple hypotheses testing (Rosetta Elucidator v3.3). The input 
for this test was the protein-level data for unenriched analysis 
(intensity for all peptides summed per sample), or the peptide-level 
data for phosphopeptides, and the raw intensities were scaled to a 
normal distribution using the Error Model in Elucidator software 
prior to ANOVA. Fold-changes and p-values are shown for all 
proteins (Table SI) and phosphopeptides (Table S2). Statistical 
cutoffs for fold-change were established for proteins based on a 
power calculation using the average biological variation within 
each group. Using the protein average %CV (23%) and 4 reps per 
group at a 95% confidence, minimum cutoffs were set to 2-fold 
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(98% powering) for proteins. ANOVA p-value of 0.05 or less was 
required. 

Western Blot Analysis 

Cryopreserved heart tissue was mechanically disrupted by 
mortar and pestie in hquid nitrogen and suspended in lysis buffer 
(1% IGEPAL CA-630, 0.5% Deoxycholate, 2% SDS, 5 mM 
EDTA in IX PBS) with protease and phosphatase inhibitor 
cocktail tablets (Roche Diagnostics, Indianapolis, IN). Samples 
were then pulse homogenized on ice via mechanical homogeni- 
zation. 

Cardiac tissue homogenates were subjected to Bicinchoninic 
acid (BCA) assay (Pierce Biotechnology/Thermo Fisher Scientific, 
Rockford, IL) for protein quantification. Western Blot for protein 
immunodetection was performed as described previously [11]. 

Metabolomics 

Organic acids were quantified using methods described 
previously [12] employing Trace Ultra GC coupled to ISQ, MS 
operating under Xcalibur 2.2 (Thermo Fisher Scientific, Austin, 
TX). 

Blolnformatlcs Analysis of Protein Data Sets 

Ingenuity Pathway Analysis (IPA, Winter 2012 Release) 
(Ingenuity Systems, Redwood City, CA) was used to classify the 
proteins according to primary function as well as for pathway 
analysis. KinasePhose 2.0 (http:/ /kinasephos2.mbc.nctu.edu.tw/ 
index.html) [13] was used to identify protein-kinase specific 
phosphorylation sites among the differentially phosphorylated 
proteins with a specificity threshold of at least 80%. 

Supporting IVlaterials and IVlethods 

Detailed sample preparation, LC-MS processing, and western 
blot methodology is described in Methods SI. 

Results 

Patient Samples 

Left ventricular tissue from 12 male patients matched for age 
and race was used in this study, with 4 patients in each of the 
following groups: ischemic failing (IF); noii-ischemic failing (NIF); 
or non-failing (NF) hearts not used for transplantation. Groups 
were closely matched for cardiac function and treatment history 
(Table 1). All patients with HF exhibited significantly lower 
ejection fractions compared to NF controls, which had normal left 
ventricular function. All HF patients received intravenous 
inotropic agents and intra-aortic balloon pump support, whereas 
only patients from the IF subset received prior coronary artery 
bypass surgery. In the NIF group, two patients were diagnosed 
with non-ischemic cardiomyopathy of unknown etiology, one 
patient had a viral cardiomyopathy, and the fourth patient 
developed HF secondary to valvular disease. 

Quantitative Analysis of Proteins in Human Left 
Ventricles 

Following the workflow in Figure lA, all heart tissue samples 
were subjected to quantitative analyses of the phospho-enriched 
and unenriched proteome. The reproducibility of the unenriched 
analytical approach was validated with a spike-in of the internal 
standard yeast ADHl digest at a known concentration. Consistent 
ADHl abundance (7% coelficient of variation (CV)) across all 12 
samples was observed (Figure IB). 
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Figure 1 . Experimental Approach. A: Illustration of the sample preparation and MS work flow. B: Reproducibility/lnternal Control. Relative levels 
of the ADHl protein spiked into each heart lysate as determined by chromatographic peak intensity. C: Venn Diagram demonstrating the protein and 
phosphoprotein yield and overlap. 
doi:1 0.1 371/journal.pone.01 041 57.g001 



Peptide/protein identification and quantification by Rosetta 
Elucidator with Mascot and IdentityE search algoritlims yielded 
expression data for a total of 850 proteins: 4,436 peptide 
annotations representing 450 proteins in the unenriched samples, 
for which relative quantitation between samples was performed 
(Figure IC; Table SI for the complete list), and site-specific 
quantitation for 823 phosphopeptides corresponding to 400 
phosphorylated proteins in the phospho-enriched samples were 
determined (Figure IC; Table S2 for the complete list). Between 
the observed unenriched proteome and the phosphoproteome, 68 
proteins overlapped (Figure IC). These 68 proteins and their 
respective fold changes and p-values are fisted in Table S3. 

Differential Expression Analysis 

Principal components analysis (PCA) was performed in order to 
observe any high-level dilferences between sample groups and to 
screen for oudier samples using the unenriched and phospho- 
enriched proteome expression data. Expression data was z-score 
transformed across all samples at the protein-level for unenriched 
samples (Figure 2A) or at the peptide level for phospho-enriched 
samples (Figure 2B), and Rosetta Elucidator v3.3 was used to 
calculate principal components. The two most prominent com- 
ponents (PGl versus PC2) are plotted in Figures 2A and 2B, and 
no significant outliers are observed for either unenriched or 
phosphoproteomes. Interestingly, in unenriched proteomes (Fig- 
ure 2A), the four NF samples separate from the failing samples 
along PC 1 ; however, in the phosphoproteome, at least three of the 
four NIF appear unique along PCI compared to the IF or NF 
groups (Figure 2B). This global observation of group classification 
suggests that the unenriched proteome may best separate faiUng 
from NF hearts, whereas phosphorylation status may play a larger 
role in distinguishing IF versus NIF. 



Using the abundance levels of significantly different proteins 
and phosphopeptides, an unbiased 2D hierarchical clustering 
analysis was performed at the protein-level (Figure 2C) or 
phosphopeptide level (Figure 2D) in order to observe the overall 
expression pattern of these molecules in each individual sample 
and to determine how these candidates segregate the failing and 
NF hearts. Figure 2C shows clear differentiation between failing 
and NF hearts, but none between the two failing etiologies was 
observed. In contrast, the phosphopeptide analyses (Figure 2D) 
show a more robust differentiation between all three groups. 

Comparison of the unenriched NIF and IF proteomes revealed 
no differential protein expression, except for carbonic anhydrase 3. 
The data from the unenriched proteome was filtered to exclude 
blood proteins and antibodies from consideration. 18 proteins 
represented by at least two high-confidence peptides and a 
significant (ANOVA, p< 0.05) fold change of at least two were 
identified between unenriched failing and NF proteomes (Table 2 
and 3). Three proteins in the IF (Table 2) and 15 proteins in NIF 
(Table 3) differed in abundance from the NF. Moreover, four 
proteins were decreased in abundance, while the remaining 
proteins were increased in abundance in the failing hearts (Table 2 
and 3). Among the proteins whose abundance was significandy 
changed in the failing human myocardial tissue were secreted 
glycoproteins, ceruloplasmin, carbonic anhydrase 1 and 3, serum 
amyloid A, and extracellular matrix proteins. Peptides from the 
fibulin family of extracellular matrix proteins (fibulin 1, fibulin2, 
and latent transforming growth factor beta binding protein 2 
(fibulin 3)) were consistently up-regulated across the NIF samples. 

Our primary goal was to determine etiology-specific changes in 
site-specific phosphorylation for the phosphopeptides. A phospho- 
peptide was considered of interest if it fulfilled the following 
criteria: 1) a ^2-fold alteration in phosphorylation state (ANOVA 
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Figure 2. Principal Components Analysis and Hierarchical Clustering. The expression data for all peptides from the unenriched (A) and 
phosphopeptide enriched (B) samples were used to perform Principal Components Analysis (PCA) after z-score transformation of the peptide 
intensities. The top two principal components are plotted in each figure, showing no extreme outlier samples among the twelve individual patients 
tested, either among the unenriched samples or phosphopeptide-enriched. The statistically-significant differentially expressed peptides for each 
experiment were used to calculate 2D hierarchical clusters in order to view sample-to-sample relationships within these differentially expressed 
signals, at the unenriched proteome level from Tables 2 and 3 (C) and from the phosphoproteome in Table 4 (D). 
doi:1 0.1 371 /journal.pone.01 041 57.g002 



p-valueS0.05) between NIF and IF; 2) significantly difTerent 
phosphorylation levels between at least one etiology of failing 
hearts and NF hearts; and 3) phosphorylation level differences 
observed were unlikely to be due to protein level differences. The 
last criteria was fulfilled by using either the unenriched protein 
levels (if available; Table S3 for the 68 overlapping proteins) or the 
comparative levels of other phosphopeptides from the protein of 
interest (Table S4) and demonstrating that the levels of these other 
phosphopeptides are not significantly different between the three 
heart tissue groups nor do they trend in the same direction as the 
phosphopeptide under consideration. 



From this analysis, 26 proteins were identified to contain at least 
one differently phosphorylated site between NIF and IF tissue that 
met our criteria (Table 4). Some of these proteins, such as lyric 
and leiomodin 1, were differentially phosphorylated in the two 
types of heart failure. Other phosphoproteins demonstrated either 
a significant hyper-phosphorylation or dephosphorylation in one 
or more amino acid sites, while phosphorylation levels at the 
corresponding amino acid sites were unchanged in the other 
etiology compared to NF tissue. 
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Table 2. Proteins with statistical significant differential expression between IF and NF heart tissue. 





Entry Name 


Protein Name 


Peptide Count Fold Change 


p-value (ANOVA) 


Protein Function 


CAH1_HUMAN 


Carbonic anhydrase 1 


2 7.24 


1.6x10"' 


Acid-Base balance 


A2MG_HUMAN 


Alpha-2-macroglobulin 


33 2.31 


1.0x10"" 


Protease Inhibitor 


SAA_HUMAN 


Serum amyloid A protein 


2 -7.32 


1.0x10"" 


Acute Phase Protein 


Data was filtered to eliminate blood proteins and antibodies and to show only proteins with Protein Prophet probability >0.8, FDR-corrected p-value <0.05, absolute 



fold-change >2. Proteins in bold were validated by western blot analysis. 
doi:l 0.1 371 /journal.pone.Ol 041 57.t002 



Validation Western Blot Analysis 

The unenriched proteomic mass spec data were validated by 
Western Blot analysis of independently prepared tissue homoge- 
nates from the same 12 patients (Figure 3). Western blot 
quantitative data is shown in Figure 4. These Western blots 
confirmed the altered changes in IF and NIF compared to NF that 
were revealed by the mass spec data: an increase in carbonic 
anhydrase, ceruloplasmin, Fibulinl, Fibulin2 and Alpha 2-HS 
glycoprotein (FETUA); a decrease in serum amyloid A; and an 
increase in 0(2-Macroglobulin (a2M). Interestingly, a cleaved 
version of a2M was the prevalent form of a2M in both NIF and IF 
hearts, whereas the full-length a2M was the prevalent species in 
NF heart tissue (Figures 3 and 4). 

Metaboiomic analysis 

Pyruvate dehydrogenase El component subunit alpha (ODPA 
or PDH) catalyzes the overall conversion of pyruvate to acetyl- 
CoA and CO2, and provides the primary link between glycolysis 
and the tricarboxylic acid cycle. A significant decrease in 



phosphorylation of PDH at Ser 300 was observed in IF heart 
compared to both NIF and NF samples (Table 4) suggestive that 
PDH is more active in IF hearts than NIF or NF hearts. Mass 
spectrometry methods were used to examine pyruvate content in 
the same 12 human myocardial samples (Figure 5). The IF tissue 
had the lowest levels of pyruvate, supportive of our hypothesis that 
of the three groups of patients, PDH is most active in IF. Another 
enzyme that influences pyruvate levels in the myocardium is 
lactate dehydrogenase (LDH). LDH catalyzes the conversion of 
pyruvate to lactate and back. Lactate levels were the highest in IF 
(Figure 5) suggestive that in IF tissues both PDH and LDH are 
more active than in NIF. 

Bioinformatics Analysis 

Focusing on differential protein or phosphoprotein levels 
between NIF and NF samples. Ingenuity Pathway Analysis (IP A) 
was used to construct an interaction network of proteins (Figure 6). 
Many of the proteins that were significantly different (> 2-fold 
change, p<0.05) were glycoproteins, proteoglycans, and structural 



Table 3. Proteins with statistically significant differential expression between NIF and NF heart tissue. 





Entry Name 


Protein Name 


Peptide Count 


Fold Change 


p-va\ue (ANOVA) 


Protein Function 


CAH3_HUMAN 


Carbonic anhydrase 3 


2 


7.02 


1.2x10"°^ 


Acid-Base balance 


LTBP2_HUMAN 


Latent-transforming growth 
factor a-binding protein 2 


6 


4.12 


9.8x10"°' 


Elastic Fiber Structure 


CAH1_HUMAN 


Carbonic anhydrase 1 


2 


4.06 


3.1 xlO""" 


Acid-Base balance 


ASPN_HUMAN 


Asporin 


9 


3.03 


1.5x10"°^ 


Cartilage Homeostasis 


C06A3_HUMAN 


Collagen alpha-3(VI) chain 


6 


2.93 


5.6x10"°' 


ECM Fibrillar Protein 


FBLN2_HUMAN 


Flbulln-2 


4 


2.52 


7.1 xlO"°' 


ECM Remodeling Protein 


MFAP4_HUMAN 


Microfibril-associated 
glycoprotein 4 


5 


2.50 


1.4x10"°^ 


Elastic Fiber Formation 


FIBG_HUMAN 


Fibrinogen gamma chain 


13 


2.29 


1.4x10°^ 


Primary Platelet Receptor 
Binding Site 


FBLN3_HUMAN 


EGF-containing fibulin-like 
extracellular matrix protein 1 


2 


2.25 


1.7x10"°^ 


Elastic Fiber Formation 


FBLN1_HUMAN 


Flbulln-1 


2 


2.16 


4.7x10"°^ 


ECM Organization 


CERU_HUMAN 


Ceruloplasmin 


15 


2.09 


4.0x10"°' 


Copper Homeostasis 


DERM_HUMAN 


Dermatopontin 


5 


2.08 


2.8x10"°^ 


Fibroblast Cell Adhesion 


MYH2_HUMAN 


IVlyosin-2 


2 


-2.08 


4.5x10"°' 


Skeletal Muscle Contraction 


ATPD_HUMAN 


ATP synthase subunit delta, 
mitochondrial 


2 


-2.22 


6.2x10"°" 


ATP Synthase Core Subunit 


SAA_HUMAN 


Serum amyloid A protein 


2 


-6.41 


1.7x10"°" 


Acute Phase Protein 



Data was filtered to eliminate blood proteins and antibodies and to show only proteins with Protein Prophet probability >0.8, FDR-corrected p-value <0.05 and 
absolute fold-change >2. Proteins in bold were validated by western blot analysis. 
doi:l 0.1 371/journal.pone.Ol 041 57.t003 
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Figure 3. Western blot validations of selected statistically significant proteins identified by proteomics analysis. Protein extracts from 
IF, NIF, or NF samples (40 |i.g) were subjected to polyacrylamlde gel electrophoresis. Replica nitrocellulose blots were incubated with anti- carbonic 
anhydrase, -ceruloplasmin, -fibulin 1, -fibulin 2, -serum amyloid A, -fetuin A, -alpha 2 macroglobulin, -PDH Serine 300, or -sarcomeric actin. 
dor:1 0.1 371/journal.pone.01 041 57.g003 



proteins, all components of the extracellular matrix (ECM). More 
ECM proteins were discovered when the inclusion criteria were 
expanded to significantly difierent proteins between NIF and NF 
with <2-fold change (min fold change 1.46) using the IPA analysis 
(Figure 6, light green). Pathway analysis also suggested that Akt, 
SMAD3, MMP14, and aryl hydrocarbon receptor (AHR), all of 
which have been implicated in cardiac remodeling [14—17], may 
be central regulators controlling the expression levels of many of 
the proteins included in the IPA analysis. IPA analysis validation 
was performed by Western blot analysis, which established that 
several of these more central proteins (AKT, SMAD3, and AHR) 
were also differentially expressed in NIF (Figures 7 and 8). 

Further bioinformatics analysis of the phosphorylation data 
demonstrated that nine of the differentially phosphorylated 
proteins are possible targets of casein kinase 2 (Ck2, Figure 9). 
The four differentially phosphorylated proteins in NIF (LYRIC, 
SHRC, FHL2, K.-XPO and MPRI) were all dephosphorylated in 
this disease state, while four of the six dififerentially phosphorylated 
proteins in IF (LYRIC, LARP7, HSP90A, ACINU) were 
increasingly phosphorylated. 



Discussion 

Heart failure remains a progressive disease with limited 
treatment options. Current standard of care does not differentiate 
between ischemic and non-ischemic etiologies, treating them 
similarly. However, as this study suggests, there are important 
protein phosphorylation differences between the two IIF (Etiolo- 
gies. Given the variety of ceUular processes these proteins are 
involved in, these differences may be important in the design of 
therapies for heart failure. 

Our analysis revealed the following general differences between 
the two main t)pes of heart failure: 1) Protein abundances were 
similar between both types of heart failure compared to NF 
controls, with one exception; and 2) Site-specific phosphorylation 
differences on many proteins clearly distinguished non-ischemic 
and ischemic HF. 

Protein abundance differences between non-ischemic 
and ischemic heart tissue 

The only protein that displayed statistically significant differen- 
tial expression between NIF and IF was carbonic anhydrase 3 
(CAH3). CAH3 levels were 2.83-fold higher in the NIF than in the 
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Figure 4. Quantitation of carbonic anhydrase 1, ceruloplasmin, fibulin 1, fibulin 2, serum amyloid A, fetuin A, and total and 
cleavage products of alpha 2 macroglobulin as determined by densitometry analysis and expressed as relative intensity (arbitrary 
units (AU)) normalized to actin levels. *: significantly elevated relative to NF samples, f: significantly elevated relative to both NIF and IF. 
doi:1 0.1 371 /journal.pone.01 041 57.g004 



IF samples. Why the levels of CAH3 increase in HF, especially in 
NIF, is unclear. In general, many tissues contain carbonic 
anhydrases, which interconvert carbon dioxide and bicarbonate 
to maintain acid-base balance and help transport carbon dioxide 
out of tissues, especially metabolically active tissues like heart and 
brain. Therefore, the increase of CAH3 could be a sign that these 
hearts are increasingly metabolically active in an attempt to 
compensate for the decrease in function. 

Protein abundance differences between failing and non- 
failing human tissue 

Extracellular matrix remodeling in the context of 
NIF. The involvement of ECM remodeling in NIF was our 
most robust fmding from evaluation of protein level differences 
between failing and non-failing tissues (Figure 6). This finding is 
not only consistent with other studies that have used proteomic 
analysis to investigate protein expression variation in cardiac 
disease, but it complements their findings as well [18,19]. The 
Barallobre-Barreiro et al study, which focused on the ECM 
remodeling in ischemia/reperfusion injury using a pig model, 
found an interaction network centered on transforming growth 
factor-pi (TGFBl), a known profibrotic molecule[19]. TGFBl 
was also heavily involved in our network (Figure 6). Additionally, 



our network demonstrated what the Barallobre-Barreiro et al 
study suggested, that fibulin 1, 2, and 3 (components of elastic 
fibers) were significantly up regulated in NIF and trended toward 
up regulation in IF. These revelations were confirmed by western 
blot analysis of fibulin 1 and 2 (Figures 3 and 4). Up regulation of 
fibulin- 1 is associated with cardiac disease, [20-22] while up 
regulation of fibulin-2 and fibulin-3, although not implicated in 
cardiac disease, is associated with a variety of other fibrotic 
disorders [23]. Down regulating fibulins may be a therapeutic 
target in future treatments of heart failure. 

Oxidative stress, cellular damage, and 
inflammation. NIF and IF samples had similarly modulated 
proteins involved in the oxidative stress and inflammatory 
pathways, including ceruloplasmin, heat shock protein 90, serum 
amyloid A, and a2-macroglobulin (oi2M). These proteins had 
comparative abundances in human heart tissue that were 
confirmed by Western blot analyses (Figures 3 and 4). 

When exposed to activated proteases, a2M exposes a "bait 
region" within its structure that traps the protease. Following 
protease binding, a2M promotes clearance and degradation of the 
bound protease while undergoing limited proteolysis [24]. This 
limited proteolysis, however, does result in a2M cleavage, and our 
findings reveal an increase in a2M cleavage products. Therefore, 
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Figure 5. Evaluation of lactate and pyruvate levels in NF, non- 
failing; IF, ischemic failing; NIF, non-ischemic failing human 
myocardial tissues. Data is reported as fluorescent intensity relative 
to NF tissue. 
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the increased levels of cleaved a2M in HF tissue could be 
protective attempts to limit the effects of harmful proteases. In 
addition to binding proteases, the native form of a2M can bind to 
damaged proteins in order to limit cellular damage and abnormal 
protein deposition [25]. However this broader protective function 
of a2M is lost upon exposure to proteases [25] and this decreased 
ability of ot2M to clear damaging protein deposits could be 
contributing to HF progression. Whether a cause of or response to 



cellular damage, the cleaved form of a-2M could be a useful 
biomarker for the progression of HF. 

NIF and IF etiologies have different phosphorylation 
signatures 

While the unenriched data set demonstrates the common 
pathways underlying both types of heart failure, the phosphory- 
lation analysis indicates differences between the two disease 
etiologies. These chflFerences can be grouped as: 1) etiology-specific 
differences in phosphorylation at a particular site; and 2) 
significantiy different single- or multiple-site phosphorylation in 
only one type of HF (compared to NF control). Another notable 
finding was that over one-third of the differentially phosphorylated 
protein sequences (9/26) were potential targets of casein kinase 2 
(Ck2) (Figure 9). Common functional themes of the significant 
phosphoproteins include cell growth/death, cardiac muscle 
development and/or function, and stress response, suggesting 
etiology-specific regulation of these functions contribute to the 
development or progression of heart failure. 

Etiology-specific differences in phosphorylation. Lyric, 
also known as metadherin, is a transcription co-factor. Increased 
Lyric content and phosphorylation play a major role in cancer cell 
survival [26,27]. Our finding is the first report indicating an 
association between phosphorylation of Lyric and heart failure. 
The NIF heart samples in this study displayed increased 
phosphorylation at lyric amino acid 298, whereas IF samples 
displayed a decrease in lyric phosphorylation at the same site 
(Table 4). Although no studies have established the role of 
phosphorylation at site 298 in the regulation of Lyric activity, 
other regions on the Lyric protein are involved in protein/ protein 
interactions and appear to be regulated by phosphorylation events. 




Figure 6. interaction network of non-ischemic failing hearts. Relevant interactions of the differentially expressed proteins and their 
relationships with selected pathologies are depicted. 
doi:1 0.1 371 /journal.pone.01 041 57.g006 
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Figure 7. Western blot analysis of central proteins identified from pathway analysis. Protein extracts from IF, NIF, or NF samples (40 \ig) 
were subjected to polyacrylamide gel electrophoresis. Replica nitrocellulose blots were incubated with anti- SMAD3, -AKT, -MMP14, -AHR or - 
sarcomeric actin. 

doi:1 0.1 371/journal.pone.01 041 57.g007 



Assuming that this phosphorylation/dephosphorylation of Lyric at 
298 also functions as a toggle switch for protein binding, the 
increased phosphorylation in IF could be an attempt to recruit a 
partner transcription factor to activate cell growth in the viable 
myocardium to compensate for the injury already suffered. 
Conversely, dephosphorylation of the pro-growth Lyric protein 
could be contributing to progression of NIF heart failure by 
activating pro-cell death pathways or suppressing pro-cell growth 
pathways. 

Leiomodin-1 is primarily expressed in smooth muscle, with 
some expression in the heart, and is speculated to play a role in 
cardiac hypertrophy and muscle contraction [28]. Phosphoryla- 
tion of leiomodin-1 was suggested to be a mediator of the actin 
cytoskeleton formation in human trabecular meshwork cells [29]. 
Multiple amino acid sites on leiomodin-1 exhibited greater 
phosphorylation in the NIF samples than the IF samples, with 
differences ranging from 4 to 9 fold greater, depending on the site 
(Table 4). The presence of leiomodin-1 in the smooth muscle and 
the cytoskeleton Knks well to a pathological role in heart failure. 
Generation of mutants in leiomodin-1 unable to be phosphory- 
lated at the amino acids revealed herein wUl be critical in 
determining the involvement of this relatively uncharacterized 
protein in the development of these two distinct types of heart 
failure. 



Significantly different phosphorylation only seen in one 
HF etiology. Phosphorylation of basigin at serine 368 near the 
C-terminus was increased in IF (compared to NIF and NF). 
Basigin is a transmembrane glycoprotein that induces MMP 
activity and subsequent ECM breakdown that is associated with 
adverse tissue remodeling seen in many pathological conditions 
such as HF[30]. Basigin is also connected to energy metabolism by 
directly mediating monocarboxylate transporter (MCT) localiza- 
tion and activity in the sarcolemma. Dissociation of basigin from 
MCTl, the primary isoform in heart, reduces cellular lactate 
influx limiting lactate oxidation inside the cardiomyocyte [31]. 
Considering basigin-MCT 1 interactions, the phosphorylation state 
of basigin could affect lactate influx into the cardiomyocyte and 
thus influence cardiac metabolism. In fact, the C-terminus of 
basigin is necessary for dimerization and association with MCTl 
[32], and a serine-to-alanine mutation at serine 193 disrupted 
basigin dimerization and MMP2 induction [33]. This data 
suggests that phosphorylation of basigin may alter its activity 
and function, and could contribute to the increase in lactate pools 
seen in the IF. The role of differential phosphorylation states of 
basigin in heart failure deserves further attention. 

POPDl, a member of the Popeye domain containing family, is 
a membrane protein abundantly expressed in the heart. POPDl 
has been shown to play an important role in cardiac pacemaker 
cells, cellular adhesion, and cellular signaling, especially in times of 
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Figure 8. Quantitation of SMAD3, AKT, MMP14, and AHR as determined by densitometry analysis and expressed as relative 
intensity (arbitrary units (AU)) normalized to actin levels. *: significantly elevated relative to non-failing, ■f: significantly elevated relative to 
both NIF and IF. 
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stress [34,35]. Some literature suggest that POPDl content is 
decreased in human heart failure [36]. While the ()\cTall 
abundance of POPD 1 was not changed between the three groups 
in our data, POPD 1 phosphorylation at 2 sites was diminished in 
NIF heart relative to NF heart, while phosphorylation at these sites 
was minimally changed in IF (Table 4). POPDl phosphorylation 
has not been examined, but this protein is important in the RhoA 
signaling pathway, a key regulator of myosin light chain 
phosphorylation [37]. Myosin light chain phosphorylation is 
decreased when POPDl is overexpressed, and a recent study 
demonstrated that the transition from compensated to decompen- 
sated heart failure is centrally mediated by a reduction in myosin 
light chain phosphorylation [38]. Interestingly, in our study, slow 
cardiac myosin regulatory light chain 2 (MLRV) was significandy 
less phosphorylated at serine 15 in NIF compared to IF heart 
(Table 4). Together, these findings suggest etiolog\'-specific phos- 
phorylation events in the RhoA/POPD 1 /MLRV pathway may 
play a role in the development of non-ischemic heart failure. 

Multiple sites on SORBS2 exhibited increased phosphorylation 
in NIF samples compared to the other two groups (Table 4). 
SORBS2 is an adaptor protein that mediates interactions between 



structural proteins, membrane proteins, other signaling molecules, 
and actin filaments in cardiomyocytes [39]. From a mass 
spectrometry based proteomics approach, SORBS2 has been 
recently identified as a protein released from cardiac tissue 
immediately following acute myocardial infarction [40] . SORBS2 
also promotes the ubiquitination and degradation of c-Abl, an 
important regulator of the actin cytoskeleton and apoptosis, [41] 
and alterations in SORBS2 activity by dilferential phosphorylation 
could activate the apoptotic cascade contributing to the develop- 
ment of non-ischemic heart failure. SORBS2 is expressed 
exclusively within the myofibril Z-bands, which link the sarco- 
meric contractile units together, and regulates the signaling 
cascade needed for proper force production and transmission in 
these contractile units [39,42]. Another sarcomeric protein that 
plays a key role in force transmission at the Z-line, titin, was 
differentially phosphorylated in NIF (Table 4). The phosphoryla- 
tion of these two structural proteins may contribute to the 
disruption of the tighdy regulated contractile apparatus leading to 
the cardiac contractile deficiencies in NIF. 

Phosphorylation/dephosphorylation of PDH/ODPA is com- 
plex. The activity of PDH is regulated via phosphorylation of three 
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Figure 9. Casein Itinase phosphorylation targets. Nine of the differentially phosphorylated proteins shown in Table 4 are possible targets of 
casein kinase. 
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serine resides (232, 293, and 300) [43]. In vitro, phosphorylation at 
a single site is sufficient to inactivate PDH, [43] and dephosphor- 
ylation by pyruvate dehydrogenase phosphatases activate PDH 
[44]. A decrease in phosphorylation of PDH at Ser 300 is observed 
in IF, while phosphorylation at Ser 293 is unchanged or only 
modestly changed. In non-failing, well-perfused heart, fatty acids 
provide 60-90% of the energy for ATP production, with the 
remaining 10-40% derived from carbohydrate (glucose and 
lactate) oxidation [45]. In the failing heart, the preferential fuel 
source switches to glucose from fatty acids. The marked decrease 
in PDH phosphorylation in IF supports the current dogma 
regarding the switch towards glucose, at least in IF, and suggests a 
mechanism by which this is occurring, namely reduction of PDH 
inhibition. Interestingly, PDH phosphorylation was not signili- 
cantiy changed in NIF hearts, suggesting that this shift in energy 
substrate is etiology specific. The trend in lower pyruvate levels in 
the IF tissue fits with our hypothesis that PDH activity is increased 
in IF. Despite the differences in pyruvate and PDH Ser 300 levels 
between HF etiologies, PDH enzyme activity was not statistically 
dilferent between groups (data not shown, p = 0.866). In all, the 
data suggest that PDH activity may differentiate between NIF and 
IF, but additional studies will be required. Whether the increase in 
glucose oxidation within the ischemic heart represents a protective 
response to chronic ischemia or a maladaptive response that 
further stresses the ischemic heart remains unclear. 

Role of Casein Kinase 2 in the development of heart 
failure: Etiology specific? 

Of the 26 differendy phosphorylated proteins, nine had sites 
that were predicted targets of casein kinase 2 (Ck2, Figure 9). Ck2 
is a ubiquitous serine/threonine protein kinases thought to have 
over 300 protein substrates that modulate a variety of cellular 
processes, including cell cycle control, cellular differentiation, and 
proliferation [46]. Classified as a "messenger-independent kinase," 



Ck2 is a constitutively active kinase that is regulated through 
protein-protein interactions, localization, and extent of phosphor- 
ylation and oligomerization [46,47]. Ck2 activity has been shown 
to both induce hypertrophy and inhibit apoptosis within cardio- 
myocytes, primarily through its interaction with apoptosis 
repressor with caspase recruitment domain (ARC), histone 
deacetylases, and the tumor suppressor gene p27 [48-52]. Many 
studies have demonstrated that Ck2 activity is critical for ischemic 
and stress-induced cardiac hypertrophic growth [48-54]. 

The majority of significant Ck2-specific sites in the IF hearts 
were increasingly phosphorylated, while all of the significant 
phosphoproteins with Ck2-specific sites in the NIF group were all 
relatively dephosphorylated. Furthermore, Lyric, a protein show- 
ing etiology-specific chflFerences in phosphorylation, is predicted to 
be a target of Ck2 [13]. 

The observation that many of these phosphorylation sites may 
be targets of Ck2 suggest that variation in Ck2 activity may point 
to a branching point between these two disease states. In ischemic 
heart failure, the initial injury and cardiomyocyte loss triggers pro- 
hypertrophic and anti-apoptotic pathways, of which Ck2 is a 
central player. This sustained cardiac hypertrophy, due in part to 
increased Ck2 activity, leads to maladaptive cardiac remodeling 
and eventually heart failure. Conversely, in non-ischemic heart 
disease, Ck2 activity is decreased, leading to increased apoptotic 
activity and suppression of pro-growth pathways. As such, 
modulation of Ck2 activity could be an etiology-specific treatment, 
although further analysis into the role of Ck2 in the development 
of ischemic and non-ischemic heart failure is needed. 

Conclusions 

This is the first study using a multi-'omics-driven systems 
biology approach to obtain a more integrated and comprehensive 
molecular assessment of the different etiologies of heart failure. 
Using an unbiased assessment of the unenriched and phosphoen- 
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riched proteomes of two HF etiologies, combined with selected 
metabolomics and Western blot analysis, the previously undocu- 
mented molecular fingerprints of the different heart failure 
etiologies were explored. These molecular differences between 
heart failure etiologies have enabled us to complement pathways 
believed to be involved in the development of HF in general. As 
our human heart repository becomes more robust in terms of 
increased specimen number and cohort size, the evaluation of 
additional post translational modifications (PTMs) and the use of 
multiple reaction monitoring [55] for more sensitive and precise 
quantitative analyses of these PTMs becomes achievable. Com- 
prehensive understanding of the proteome of this common disease 
through study of human cardiac tissue from a robust biorepository 
win lead to better, more focused treatments, and ultimately 
improved patient outcomes. 

Supporting Information 

Excel File SI The following information is contained 
within Excel File SI. Key sheet-Contains an index of the 
remainder of the excel file sheets, sample identifiers, patient ID, 
disease etiology, peptide modifications, and Modloc scoring 
criteria. ST 1 -peptide expression data from unbiased proteomics. 
ST2-protern expression data from unbiased proteomics. ST3- 
Phosphopeptide expression data from phospho-enriched proteo- 
mics. ST4-DifiFerential Unbiased Protein Expression (Fold Chang- 
es, ANOVA) for unbiased proteomics. STS-Differential Phospho- 
peptide Expression (Fold Changes, ANOVA). 
PCLSX) 
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